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Abstract 

We have shown that the dynamics of the scalar field </>(x) = Ll G~ 1 {x) n in Brans- 
Dicke theories of gravity makes the surface area of the black hole horizon oscillatory dur- 
ing its dynamical evolution. It explicitly explains why the area theorem does not hold 
in Brans-Dicke theory. However, we show that there exists a certain non-decreasing 
quantity defined on the event horizon which is proportional to the black hole entropy 
for the case of stationary solutions in Brans-Dicke theory. Some numerical simulations 
have been demonstrated for Oppenheimer-Snyder collapse in Brans-Dicke theory. 
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1 Introduction 



Scheel, Shapiro and Teukolskyjl]] have numerically shown that Oppenheimer-Snyder collapse 
in Brans-Dicke theories of gravity produces black holes that are identical to those of general 
relativity in the final stationary stage, but behave quite differently during dynamical evolu- 
tion. For example, in general relativity the apparent horizon of a black hole is always inside 
the event horizon and the total surface area of the event horizon never decreases in any 
classical process provided that null energy condition for matter fields(i.e., T a t,k a k b > for all 
null vectors k a ) and cosmic censorship conjecture are satisfied^. In Brans-Dicke theories of 
gravity, however, they found that there are some initial epochs in Oppenheimer-Snyder col- 
lapse during which not only the apparent horizon passes outside the event horizon, but also 
the surface area of the event horizon decreases in time. Thus, the area increase theorem for 
black holes in general relativity does not hold in Brans-Dicke theory. As mentioned briefly 
in Ref. [|l[], these different behaviors are possible because the null convergence condition, i.e., 
R a bk a k b > for all null vectors k a , is violated for a dynamical black hole in Brans-Dicke 
theory. At the present paper, we investigate how the behavior of the auxilliary scalar field 
4>(x) = u G~ 1 (x)" in Brans-Dicke theory causes this violation and the oscillatory behavior of 
the surface area in detail. 

One may think that the violation of the area theorem for black holes in Brans-Dicke (BD) 
theory causes some problem to the black hole thermodynamics in BD gravity since the black 
hole entropy in Einstein gravity is propotional to the surface area of the horizon and so the 
area theorem automatically serves as the classical second law in black hole mechanics. We 
show that there exists a certain quantity defined on the event horizon which never decreases 
during any classical process in BD gravity. Moreover, this quantity coincides with the entropy 
in the case of stationary black holes in BD theory. Thus, the black hole thermodynamic 
second law in BD theory is established. Since BD theories of gravity are dynamically related 
to other theories of gravity such as higher curvature theories of gravity and dilaton gravity, 
the results for black holes in BD theory also give some insights for those in other theories of 
gravity. 

In sec. 0, we explain why black holes in BD theory behave differently from those in 
Einstein gravity by investigating the behavior of the auxilliary scalar field <p(x). In sec. |3], a 
quantity is constructed on the event horizon from various points of view and proved that it is 
always non-decreasing for arbitary dynamical processes. Finally, the validity of assumptions 
on the scalar field (f)(x) used in the proof is discussed. General feature of the proof and 
possible application to other theorems are also mentioned briefly. 

2 Black hole area in Brans-Dicke theory 

Let us consider the change of the total surface area of the horizon of any black hole along 
the null congruence of the horizon generators orthogonal to the spacelike cross-section H: 




(1) 
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Figure 1: Behavior of the expansion 9 in general relativity. 



where A is the affine parameter of null geodesies whose tangents are k a , i.e., fc a V a = d/dX, 
and 9 = d(\n\/h)/dX = V ' a k a is the expansion of the horizon generators. One can easily 
see that, if 0(A) > everywhere on the horizon and at any point in A, dA/dX > and so 
the surface area is non-decreasing always. Otherwise, it either still increases or decreases 
depending on the value of integration of 9 over the whole event horizon 0. 

Now let us consider the evolution of expansion along a hypersurface orthogonal null 
congruence of geodesies in general. It is determined by Raychaudhuri equation 

^ = -\ e2 -° 2 - R ^ kakh ( 2 ) 

where a 2 is the square of the shear. In the case that R a bk a k b > for any null vector k a , 
d9/dX is negative-definite and 

Or 

— > -■ (4) 
dX ~ 2 K J 

Suppose one has 9 < at some point on a null geodesic. Then the above inequality equation 
tells that 9~ l should meet zero within some finite affine parameter. That is to say, 9 reaches 
— oo at some finite affine parameter, giving a conjugate point as can be seen in Fig. [I]. 

When applied to black holes in a strongly asymptotically predictable spacetime(thus, no 
naked singularity exists outside the black hole region), this behavior of expansion leads that 
9 cannot be negative, i.e., 9 > everywhere on the horizon and at any point in A[@; a similar 
derivation will be explained more detaily in sec. |3|. It, therefore, shows dA/dX > proving 
the area theorem: the total surface area never decreases in any classical process. 

For black holes in Einstein gravity, the null convergence condition is equivalent to the 
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Figure 2: Oscillatory behavior of the scalar field 4>(x) in Brans-Dicke theory. 



null energy condition since field equations are 

Rab - ^Rg a b = SnGT ab . (5) 
For black holes in BD gravity, however, R ab k a k b becomes indefinite as follows; 

£ B D = - ^V a 0V a 0) + £ ma tter. (6) 

167rG 



Field equations are 

1 1 

(R ab - -Rg ab )<P = 87iGT ab + ^(V a 0V 6 - -g ab V C 0V C 0) + V a V fe - g ab V c V c <p. (7) 
And so 

R ab k a k b = 8irG(fi~ 1 T ab k a k b + cj0" 2 (rv a 0) 2 + (j)~ l k a k b V a V 6 </>. (8) 

One sees that, even if null energy condition and positivity of and u are satisfied, R ab k a k b 
is indefinite since the last term, (f)~ 1 k a k b V ' a Vb4> = 4>~ 1 k a 'Va(k b 'V b (J)) = 4>~ l d 2 (j) / d\ 2 = </>~ V; 
could be strongly negative. Now Eq. (§) becomes 

9' = - \-9 2 + a 2 + 8irG(l)- 1 T ab k a k b + cu(0'/0) 2 ] - <jr (9) 

Therefore, if the auxilliary scalar field <p(x) behaves as in Fig. ^| during the dynamical evo- 
lution of black holes, the expansion 6(X) possibly behaves as in Fig. |] along the horizon 
generators. In other words, before 9 hits the negative infinity after Ai in Fig. ||, the neg- 
ativeness of 0" starts to increase 9 up over zero. Consequently, the total surface area of 
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Figure 3: Possible behavior of the expansion 6 in Brans-Dicke theory corresponding to that 
of 4> in Fig. 0. 



the horizon is increasing until A = Ai, decreasing between Ai and A2, increasing again after 
A = A2, and finally approaches constant. 

For Oppenheimer-Snyder collapse in BD theory where an initially uniform particle distri- 
bution within a sphere is allowed to collapse to a final state, the numerical results obtained 
in Ref. |T| demonstrate this behavior described above. That is, Figure 10 in Ref. Q shows 
the indefiniteness of R a bk a k b for different values of u, and Figure 9 in Ref. [l|] shows that 
the surface area is strongly oscillating for small u and becomes monotonically increasing for 
large uo. A monotonic behavior for large u appears because, even though BD theory does 
not completely reduce to Einstein gravity in general as u — > 00, it does for physical cases as 
pointed out in Ref. Figure 9 in Ref. also shows that apparent horizons pass over the 
event horizons during epochs of 9 < 0. It is simply because the expansion of null geodesies is 
negative on the event horizon for those periods. Then trapped surfaces can exist outside the 
event horizon in contrary to black hole solutions in general relativity. The outer boundary 
of this trapped region where expansion vanishes will be the apparent horizon. Interestingly, 
"outgoing" null geodesies from a trapped surface outside the event horizon do not hit the 
singularity. Instead they again expand afterwhile perhaps due to the last term in Eq. (§) 
and escape to the null infinity finally. 



3 A non decreasing quantity 

In the thermodynamic analogy of classical black hole mechanics in Einstein gravity, the total 
surface area of the event horizon plays a somewhat special role. It is proportional to the 
black hole entropy, which has been later justified by considering quantum effect on matter 
fields around the black hole. Hawking's area increase theorem then establishes the classical 
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second law of black hole mechanics in Einstein gravity. In the previous section, however, we 
have seen that the area theorem does not hold for black holes in Brans-Dicke theory. Does 
this mean that the second law does not hold for black holes in BD theory? This question 
turns out to be irrelevent since in general "area-entropy" relation no longer holds in other 
theories of gravity Q. Indeed this modification of "area-entropy" relation happens for black 
holes in BD theory. 

The entropy for stationary black holes in Einstein gravity is given by 

S E = A/ AG = — <f d 2 xVh. (10) 
4G Jh 

Since Newton's constant G becomes dynamical in BD theory, a natural candidate for the 
entropy of stationary black holes in BD theory will be 



Sbd = ~^{ d 2 xVh<j)(x). (11) 
4G Jh 



Note that we defined <f>(x) to be dimensionless. This quantity is indeed the black hole 
entropy for stationary solutions in BD theory evaluated by many other methods such as 
Noether charge method, field redefinition, and Euclidean method^]. 

The formula above can also be obtained by considering a "physical process" derivation of 
the first law in black hole thermodynamics introduced by Wald@ and extended to a certain 
class of higher curvature theories in Ref. [[?]]. Consider a black hole stationary initially, being 
perturbed by small amount of matter falling, and finally settling down to a stationary state 
again. Small amount of matter falling into a stationary black hole, AM = J Q Ta^dTP and 
A J = — § T ab ip a dY!' in the asymptotically flat region, produces some change at the event 
horizon as follows 

AM - QAJ = [ T abX a di: b = - I T abX a dZ b = « f T ab k a k b XVhdXd 2 x (12) 

JO JH JH 

where \ a = £ a + Q(p a is the Killing generator of the horizon in leading order of the pertur- 
bation and H = 7i <8> R the black hole horizon; see details in Ref. . From field equations 
in (0) one obtains in leading order 

/ d 2 xdXVhXT ab k a k b = ^— { d 2 x / A/ dXVhXUR ab k a k b - -cf)' 2 - <P") 
jh 8nG Jh J\i cp 

~ / d 2 x / A/ dXVhX(-(f)6' - <p") 

87rG jh JXi 

— o~7^ / d 2 x[(Vhcf>) |Jf - VhX(<j>6 + 0') |*']. (13) 
o7rG Jh 

Since 9 = <p' = on both initial and final stationary stages, one finds 

AM - QAJ ~ d 2 x{Vh<P) |£ = ^ASbd, (14) 
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giving the formula in ( |TTD for black hole entropy in BD theory. 

Now let us assume that the quantity in (|TTJ) is the black hole entropy even at any moment 



of dynamical evolution of the horizon and see how it behaves. To be more general in the 
proof, let us consider the following quantity 



X — I d 2 xVhe p (15) 
G Jn 



S ~AGjh 

where e p is a scalar function locally defined on the horizon. As before, the change of S(X) 
along the null congruence generating the event horizon is 

^L = -Ll d 2 xVhe p ~6 (16) 
dX AG In 

with 

e = e + d xP . (17) 

Now let us apply the method developed in Ref. which basically follows Hawking's proof of 
the area theorem, with 9 in place of 9. Again the question is whether or not there can exist 
a point along the null geodesies at which 9 becomes negative. The Raychaudhuri equation 
shows 

d x 9 = -~0 2 - a 2 - u 2 - [R ah - V a V b p)k a k b . (18) 
If the last term is positive-definite, then one has 

d x 9 < - X -9 2 , (19) 

or 

dxt 1 > \{9l~9f. (20) 

Now suppose 9 < at some point on the horizon. Then in a neighborhood of that point one 
can deform a spacelike slice of the horizon slightly outward to obtain a compact spacelike 
surface £ which enters J _ (X + ) and has 9 < everywhere on E, being defined along the 
outgoing null geodesic congruence orthogonal to E. If cosmic censorship is assumed, then 
there is necessarily some null geodesic orthogonal to £ that remains on the boundary of the 
future of £ all the way out to X + 0. In other words, this geodesic has no conjugate point 
between £ and X + and is future complete. However, this is impossible for the following 
reason. Asymptotic flatness implies that p — > like A -1 at infinity for the case of ( p"T| ) 
and so does 9, where A is the affine parameter along an outgoing null geodesic. Therefore 
9/9 — » 1 + 0(A _1 ), so the inequality fl2T| ) implies that, as one follows the geodesic outwards 
from E, 9 reaches — oo within some finite affine parameter. Since 9 = 9 + d\p, this means 
that either 9 or d\p goes to — oo. In the former case we have a contradiction, as in the area 
theorem, since it implies there is a conjugate point on the geodesic. In the latter case it 
leads to a naked singularity outside the horizon, contradicting to cosmic censorship, as shall 
be shown for BD theory below. 
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Figure 4: Behavior of the surface area of the horizon (solid line) and black hole entropy 
(dotted line) in numerical simulation for Oppenheimer-Snyder collapse in Brans-Dicke theory 
with <jj — 0. Figure courtesy of Mark A. Scheel and Richard O'Shaughnessy. The figure was 
produced using the numerical code described in PRD51, 4208 (1995). 



Therefore, one finally proves the entropy increase theorem for black holes in BD theory 
by checking whether or not the positivity of (R ab — ^ dS bP)k a k h is satisfied. Since e p = 4>(x), 
one has 

(R ab - V a V b p)k a k b = R^k'-^k^VaV^+ik-V^/^ 2 

= %-KG(f)- l T ab k a k b + (l + uj)(k a V a <j)/<j)) 2 (21) 

which is manifestly non-negative provided that null energy condition for matter fields and 
positivity of <fi and u are satisfied. Since d\e p = /c a V a = e p d\p = <pd\p and is assumed 
to be positive^ 0), the divergence of d\p implies a curvature singularity from the equation 
(§) unless it is cancelled by other terms on the right hand side of the equation. Note finally 
that the proof described above does not require the existence of regular event horizon as in 
Hawking's area theorem. A numerical simulation for Oppenheimer-Snyder collapse in BD 
theory with uj = in Fig. f| demonstrates that the quantity defined in Eq. ( |i~TD is always 
increasing even if the surface area of the event horizon is oscillatory. From the behavior of 
the surface area of the horizon one sees that 9 is positive initially, zero near t = 60, becomes 
negative, and increases to zero finally. In order for 9 to increase from some negative value 
near the final stationary stage, the only possible way is that <p" is negative for that period 
and dominant in the equation @. In fact 0" turns to be negative near t = 75 in Fig. [5| For 
other values of uo one can refer to Figure 12 in Ref. |IJ. The black hole area in conformally 
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Figure 5: Positivity and oscillatory behavior of the scalar field <f>(x) corresponding to Fig. f|. 
Figure courtesy of Mark A. Scheel and Richard O'Shaughnessy. The figure was produced 
using the numerical code described in PRD51, 4208 (1995). 



related Einstein frame is indeed the same as black hole entropy in BD theory as shall be 
shown in Eq. ( p6|) below. 

4 Discussion 

In summary we have shown that the decrease of black hole area during the dynamical epoch of 
Oppenheimer-Snyder collapse in BD theory is due to the oscillatory behavior of the auxilliary 
scalar field 4>(x) = "G _1 (x)" in strong gravitational field region. We also have proved that 
the quantity defined in (|TTD , which is indeed the entropy of stationary black holes in BD 
theory, never decreases under arbitrary dynamical evolutions. Thus, it establishes a second 
law of black hole thermodynamics in BD theory. 

Results for dynamical black holes in BD theory also give some hints on the behavior of 
black holes in other theories of gravity. For example, let us consider curvature scalar squared 
theories of gravity 

J = / d 4 x^[-^—(R + aR 2 ) + CM,g)]. (22) 

J 1D7TG 

This higher curvature theory is dynamically equivalent to the following theory with no higher 



curvature term but with one auxilliary field 



h= Jd i x^[ T ^(l + 2am-j^$ 2 + £rn(ip,g)\. (23) 

By redefining the scalar field 1 + 2a<3> = 9 with q 2 = 1 + 2co>/3 and conformally scaling the 
metric g a b = <f) l ~ q g a b, one finds that I\ has the form of BD theory with a potential term and 
unconventional couplings between the scalar field and matter fields as follows 

h = J fay/^-^MR - ^V a 0V a 0) - 1/(0) + 2 M£ m (V, <?-*§)) (24) 

where V(4>) = 2 (1 — <f)~ g ) 2 /l6irGa. Further conformally transforming g a f, = _1 g a fc and 



redefining the field (f = (3 1 ln0 with (3 = yl6nG/(3 + 2uj), one reaches to the form of 
ordinary Einstein gravity || || 

h = ( d*x^pg\-^—R - ^VapVV - e- 2 ^V(<p) + e~ 2 ^C m (^, e~^g)\. (25) 
J v 107rG 2 

Now let us ask how the black hole area in curvature squared theories behaves. First of all, 
notice that the null convergence condition is not satisfied in this theory as well. Thus, one 
expects that the black hole area will be oscillatory durinng dynamical epochs as in the case 
of BD theory. The relationship between I and I2 shows it more explicitly. Assume an initial 
Cauchy surface such as that of Oppenheimer-Snyder collapse in I where spacetime is almost 
flat, e.g., R ~ 0. Since $ = R in I±, ~ 1 and V(<p) ~ on the corresponding Cauchy 
surface in J 2 . From the results for BD theory above, one then expects that the behavior of 
<p is oscillatory during dynamical epoch of the collapse. The "extended" null convergence 
condition in ([H]) to prove the increasement of Vh(p is still satisfied in I2 even in the presence 
of the potential V(0) and the unconventional couplings between <p{x) and matter fields as 
long as <p and u are positive. Thus the quantity Vh(p = y/h(p q on the event horizon never 
decreases. Then the "area" element \fh = (Vh(p)/(p q will be in general oscillatory for black 
holes in J unless the ratial increasement of Vh<p is bigger than that of (j) q for the increasing 
period. However, the positivity of <f)(x) is guaranteed only for a > 0. It follows because the 
form of V(<p) confines the field <fr i n the positive region only if a > since the potential barrier 
increases exponentially as (p — > —00 in ^(i.e., — ► in When a < 0, the potential 
falls down exponentially, easily becomes negative, and then the conformal factors become 
singular at = 0. It also has been shown that the theory Iq has the well-posed initial value 
formulation only if a > 0||T0| . 



As byproducts of the relationships above one finds 

5 BH = A = ^£/^=^£/w/w> 



= — I d 2 xVh{l + 2a$) = \zl d 2 xVh(l + 2aR). 
4G J Hi 4G Jhq 



(26) 
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Therefore, one sees that Hawking' s area theorem is transferee! to the increasement of black 
hole entropy in BD theory as well as in higher curvature theories through conformal trans- 



formations and an introduction of a scalar field [11]. The anaysis above can also be extended 
to more general class of higher curvature theories such as actions polynomial in -R0. 

In the proof given in sec. 0, we assumed the positivity of the scalar field <f>(x) and the 
coupling constant u in BD theory. The positivity of to is natural since, otherwise, it gives 
unphysical negative energy matter in the theory. However, the positivity of <j>{x) at any 
moment is highly non-trivial in BD theory and should be guaranteed in order for the quantity 
in (pLlf) to be interpreted as "entropy," which is positive by definition in statistical mechanics. 
From field equations in fl7|) we see that any constant can be a solution. Thus, each set of 
solutions for a given = const. (which could be negative as well) in BD theory with large uj 
gives Einstein-like gravity. For Einstein gravity this constant is determined by considering 
Newtonian limit. It gives a positive constant = 1 in the unit of G~ l and so the black hole 
entropy ( |TUp in Einstein gravity is always positive. In BD theory, however, (j){x) is highly 
dynamical and so the quantity in (|TTD could be negative at some stage, which is problematic 
to be an entropy of dynamical black holes in BD theory. We are in fact interested in a 
physical system in which the matter is initially distributed in an almost flat space. In other 
words, 4>(x) is initially positive near the unit. Now the question is whether or not there 
is a Cauchy surface in the future evolution where the field <f>(x) passes zero and becomes 
negative. As can be seen in field equations (0), it implies a curvature singularity unless 
the total energy-momentum tensor on the right hand side cancels this singular behavior. 
Therefore, — > corresponds to a naked singularity on the event horizon or outside the 
horizon, violating the assumption of cosmic censorship in BD theory which we used in the 
proof. To see whether the dynamics of the field can prevent this singular behavior, let us 
examine the relationship between BD theory J 2 and Einstein gravity 1%. In the theory J 3 , 
= corresponds to ip — > — 00. In general, the dynamics of the field ip is determined by the 
potentials coming from couplings with matter Lagrangian as well as V(<p). If the net effect 
of potentials gives rising up as (p — > — 00, the field tp cannot reach to — 00. One example is 
the higher curvature squared theory with a > through conformal rescalings as explained 
above. For BD theory without potential terms, that is, V — and q — 1 in J 2 , the rising 
up of potantial barrier as (p — > —00 is not guaranteed in general. For example, suppose a 
typical matter field, C m (if),g) ~ — |V a , 0V a "0 — V(ip). Then, 

e~ 2 ^C m ^, e-tog) ~ -[i(V^V»e"^ + e'^V^)}. (27) 

Since V a ?/>V a ?/> could be negative in general, the potential could fall exponentially, leading 
to (p — > —00 and so —>■ 0. As far as we know, however, we have not found any physical 
system in the literature which shows vanishing of <p(x). Numerical simulation for the case of 
Oppenheimer-Snyder collapse in Fig. [| also shows that the scalar field 0(x) is always positive 
during its evolution. 

Finally, let us extract some general feature of the proof shown in sec. [5| One of the most 
important behaviors of expansion 6 used in proving many global properties of spacetime 
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such as singularity theorem and area theorem is that 9 reaches — oo within some finite affine 
parameter if 9 is negative at some point and convergence condition(i.e., R a t,v a v b > for 
all non-spacelike vectors v a ) is satisfied as explained in sec. |2] briefly. In general relativity 
the convergence condition is satisfied if some suitable energy condition for matter fields is 
assumed. In other theories of gravity, however, this is not true any more as seen, for instance, 
in sec. |2]. Thus, it is not clear whether or not some global theorems in general relativity or 
modified forms of them where the convergent behavior of 9 is used for proof still hold in other 
theories of gravity. In this paper we have shown that a modified "expansion" 9 plays the 
same role as 9 under same energy conditions for matter fields. For example, one may define 
a modified "trapped" surface as a compact, two dimensional, smooth spacelike surface where 
all "ingoing" and "outgoing" null geodesies orthogonal to it have negative 9, instead of having 
negative 9. In fact, for BD theory, this 9 is nothing but the expansion in the conformally 
related Einstein frame with an overal multiflication factor. However, the general feature of 
the proof in this paper may still be applicable to other class of gravity theories which are not 
conformally related to Einstein gravity. In addition, it also gives some hints on investigating 
how robust many global properties of spacetime in general relativity are under the change 
of dynamics of gravitational fields. In other words, by using the convergent behavior of 9, 
we may extend the validity of many global results in general relativity to other theories of 
gravity as done for the area theorem in this paper. 

The author would like to acknowledge useful discussions with J. Samuel, T. Jacobson, J.H. 
Cho, and J.H. Yoon. Especially, I would like to thank M.A. Scheel and R. O'Shaughnessy 
for their comments and kindly making diagrams used in this paper. 
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